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Abstract 

Some tools and ideas are interchanged between random matrix theory and multi- 
variate statistics. In the context of the random matrix theory, classes of spherical 
and generalised Wishart random matrix ensemble, containing as particular cases the 
classical random matrix ensembles, are proposed. Some properties of these classes 
of ensemble are analysed. In addition, the random matrix ensemble approach is ex- 
tended and a unified theory proposed for the study of distributions for real normed 
division algebras in the context of multivariate statistics. 



1 Introduction 



Random matrices first appeared in studies of mathematical statistics in the 1930s but did 
not attract much attention at the time. However, in time the foundations were laid for one 
of the main areas of present-day statistics, Multivariate Statistics or Multiv ariate Stati stical 
Analysis. The a dvances made have been compiled in ex cell ent books such as Rovl ( 19571 ) and 
Anderson! ( 19581). l ater inlSrivastava and Khatri ( 1979f ) and lMuirheadl (|l982f) and recently in 
Gupta and Nagari (|2000h and lEatonl (|2007l ). among many other important texts. In general, 
these writings are based on the multivariate normal (Gaussian) distribution. Nevertheless, 
over the last 20 years, a new current has been developed and consolidated to constitute what 
is now known as Genaralised Multivariate Analysis, in which the multivariate normal distri- 
bution is replaced by a family of elliptical distributions, containing as particular cases the 
No rmal, Student-t, P e arson type VII, Logistic and Kotz distributions, among many others, 



see 



Fang et all (|l990l) . iFang and Zhang! (|l990h and iGupta and Vargal (|l993T ). This family 



contains many distributions with heavier or lighter tails than the Gaussian distribution, 
thus offering a more flexible framework for modelling phenomena and experiments in more 
general situations. 
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Random matrices are not only of interest to statisticians; they are studied in many disci- 
plines of science, engineering and finance. This approac h has been consol i dated and is gener- 
ally kn own a s Random Matrix Theory, see Metha ( 1991 ). Forrester] ( 20091 ), [Ratnarai ah et all 
(2005a)) and lEdelman and Raol ( 20051) . among many others. Fundamentally, this theory is 
concerned with the following question: consider a large matrix whose elements are random 
variables with given probability laws. Then, what can be said about the probabilities of 
a few of its eigenvalues or of a few of its eigenvectors? This question is relevant to our 
understanding of the statistical behaviour of slow neutron resonances in nuclear physics, a 
field in which the issue was first addressed in the 1950s and where it is intensively studied 
by physicists. The question was later found to be important in other areas of physics and 
mathematics, such as the characterisation of chaotic systems, the elastodynamic properties 
of structural materials, conductivity in disordered metals, the distribution of the values of 
the Riemann zeta function on the critical line, the enumeration of permutations with certain 
particularities, the counting of certain knots and link s, qu antum gravity, q uantum chromo 
dynamics, string theory, and others, see iMethal (|l99ll ) and iForresten (|2009h . This approach 
has also motivated the study of random linear systems (operators) , random matrix calculus 
(calculus of Jacobians of matr ix factorisations) and num erical stochastic algorithms, among 
other mathematical fields, see Edelman and Raol (2005). A very interesting characteristic of 
the results obtained in random matrix theory is that they are vaild for real, complex and 
quaternionic random matrices. 

From the standpoint of the multivariate statistics, the main focus of statistical studies 
has been developed for real random matrices. However, many texts have proposed the 
extension of some of these tests to the case of c omplex random matrices, see for example 
Woodins) (Il956l) iGoodmar] (Il963h.l,lamesl (I1964T) iKhatril (|l965l ). lHannanl (|l970l Section 6.6, 



p. 295). iMethail (|1997l ) and lMicheas et all ([2006), among many others. St udies examinin, 
the case of quaternion ic random matrices case are much less common, see iBhavsarl ( 200 
and lLi and Xuel (|2009h . 

Form the point of view of the statistical and random matrix theory, the relevance of 
the octonions is not clear. A n exce l lent r eview of the history, construction and many other 
properties of the octonions is IBaezI (|2002l) . where it is stated that: 



"Their relevance to geometry was quite obscure until 1925, when Elie Car- 
tan described 'triality' — the symmetry between vector and spinors in 8- 
dimensional Euclidian space. Their potential relevance to physics was no- 
ticed in a 1934 paper by Jordan, von Neumann and Wigner on the founda- 
tions of quantum mechanics... Work along these lines continued quite slowly 
until the 1980s, when it was realized that the octionions explain some cu- 
rious features of string theory... However, there is still no proof that 
the octonions are useful for understanding the real world. We can 
only hope that eventually this question will be settled one way or another. " 

Others problems, such as the existence o f real or imaginary e i genva lues for 2x2 and 3x3 
octonionic Hcrmitian matrices are treated in lDrav and Manogue (1998) and lDrav and Manogud 
(1999). For the sake of completeness, in the present study the case of the octonions is con- 
sidered, but the veracity of the results in this case can only be conjectured. In particular, 
Forrester (2009, Section 1.4.5, pp. 22-24) it is proveed that the bi-dimensional density func- 
tion of the eigenvalue, for a Gaussian ensemble of a 2 x 2 octonionic matrix, is obtained from 
the general joint density function of the eigenvalues for the Gaussian ensemble, assuming 
m = 2 and /3 = 8, see Section [3] 

In the present study, we propose an interchange between some ideas and tools of random 
matrix theory and multivariate statistics. Specifically, in Section [2] general expressions are 
given for the Jacobians of certain matrix transformation, together with the factorisation for 



2 



real normcd division algebras, a question that is fundamental in the study of general ran- 
dom matrix variate distributions. Vector-spherical, spherical, vector-spherical-Laguerre and 
spherical-Laguerre random matrix ensembles are studied in Section [31 where it is observed 
that the classical random matrix ensembles are obtained as particular cases of the above 
general random matrix variate distributions. Finally, the density function of the matrix 
variate normal, Wishart, elliptical, generalised Wishart and beta type I and II distributions 
are studied for real normed division algebras. 



2 Jacobians 

In this section we discuss the formulas containing the Jacobians of familiar matrix trans- 
formations and factorisations using the parameter (3 (defined below). First, however, some 
notation must be established. 

The parameter (3 has traditionally been used to count the real dimension of the underlying 
normed division algebra. In other branches of mathematics, the parameter a = 2//3 is used. 



Table 1: Values of (3 — 2/ a parameter. 






a 


Normed divison algebra 


1 


2 


real (3?) 


2 


1 


complex (£) 


4 


1/2 


quaternionic (Sj) 


8 


1/4 


octonion (£)) 



Remark 2.1. In iBaea ( 20021 ) the following remarks are made: 



"There are exactly four normed division algebras: the real numbers (3i), 
complex numbers (€), quaternions (Sj), and octonions (D). The real num- 
bers are the dependable breadwinner of the family, the complete ordered 
field we all rely on. The complex numbers are a slightly flashier but 
still respectable younger brother: not ordered, but algebraically complete. 
The quaternions, being noncommutative, are the eccentric cousin who is 
shunned at important family gatherings. But the octonions are the crazy 
old uncle nobody lets out of the attic: they are nonassociative. " 

In addition, take into account that 5ft, £, ft and O are the only normed division algebras; 
moreover, they are the only alte rnative divi sion algebras, and all division algebras have a 



real dimension of 1,2,4 or 8, see IBaea (|2002l Theorems 1, 2 and 3). 



Let £^ n be the linear space of all n x to matrices of rank m < n over J with to distinct 
positive singular values, where J denotes a real finite- dimensional normed division algebra. 
Let g" nxm be the set of all n x m matrices over The dimension of fi nxm over 5ft is (3mn. 
Let A 6 £ nxm , then A* = A denotes the usual conjugate transpose. The set of matrices 
Hi 6 ^nxm suc j 1 ^hat H^Hi = I m is a manifold denoted n , termed Stiefel manifold (Hi 
is also known as semi-orthogonal ((3—1), se mi-unitary (0 = 2), semi-s ymplectic (f3 = 4) 
and semi- exceptional type (j3 = 8) matrices, see lDrav and Manoguel (|l999f) ). The dimension 
of n over 5ft is \J3mn — m(m — l)/3/2 — to]. In particular, V^ n m with dimension over 5ft, 
[m(m + l)/3/2 — to], is the maximal compact subgroup UP(m) of m and consist of all 
matrices H G g m >< m such that H*H = I m . Therefore, il^(m) is the real orthogonal group 
O(m) ((3 = 1), the unitary group U(m) ((3 = 2), compact symplectic group Sp(m) ((3 = 4) 
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or exceptional type matrices Oo(m) ((3 — 8), for J = 5ft, C, Sj or D, respectively. Denote by 
6^ the real vector space of all S G y™ xm sucn that S = S*. Let be the cone of positive 
definite matrices S G g™ xm ; then is an open subset of 6^. Over 5ft, 6^ consist of 
symmetric matrices; over (t, Hermitian matrices; over f), quaternionic Hermitian matrices 
(also termed self-dual matrices) and over O, octonionic Hermitian matrices. Generically, 
the elements of are termed as Hermitian matrices, irrespective of the nature of 
The dimension of over 5ft is [to(to — l)/3-t-2]/2. Let be the diagonal subgroup of m 
consisting of all D G ^mxm^ -p _ dia,g(di, , . . j d m ). Let 1^(to) be the subgroup of all upper 
triangular matrices T G ^ mxm such that tij = for 1 < i < j < to; and let T^(to) be the 
opposed lower triangular subgroup T^(to) = ^T^(to)^ . For any matrix X G j« x "\ ^-jc 

denotes the matrix of differentials (dxij). Fi nally, we define t he measure or volume element 
(dX) when X G y mx ™, 6&, or V£ „, see IPimitriul (120021 ). 



If X G ^ nXm then (g?X) (the Lebesgue measure in J lixm ) denotes the exterior product 
of the fimn functionally independent variables 

n m /3 
i=l j = l fe=l 

Remark 2.2. Note that for Xij G # 

/3 



c?Xij = A dx 



(fc) 

y ■ 

fc=i 



In particular for ^ = 5ft, £, .f) or O we have 



• x^ G 5ft then 



• = xf + ixf) G £, then 



fe=i 



2 

dxy = cfecg 3 A dxf = f\ dxf 

k=l 



• Xij = xf + ixf + jxf + kxf G io, then 



4 

da^ = dxf A dsg> A da;g> A dxf = f\ dxf 



k=l 

(1) , (2) , (3) . (4) . (5) . (6) . (7) . (8) ^ 

x\j + eixlj' + e 2 x\j + e^xlj + e 4 x i y + e^x^' + e§x\j + ejx^' G £>, then 



dXij = dxf A dxf A dxf A dxf A dxf A dxf A dxf A dxg } = f\ dxf 

k=l 



If S G (or S G T^(to)) then (dS) (the Lebes gue measure in 6^ or in T^(to)) denotes 
the exterior product of the m(m+ l)/3/2 functionally independent variables (or denotes the 
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exterior product of the m(m — 1)0/2 + n functionally independent variables, if su € 3? for 
alH = 1, . . . , m) 

m ft 

A A *>$\ 

(dS) = I k=1 „ 

m rri p 

f\d Sii f\ l\ds§\ if a* eft. 

2 — 1 ?<j fc=l 

Remark 2.3. Since generally the context establishes the conditions on the elements of S, 
that is, if Sij E £ £, E f) or E O. It is considered that 

m P m m P 

( dS ) = A A d i fe) -A^A A rf i fc) - 

i<j k=l i=l i<j fe=l 

Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S G 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive matrix). In 
the re a l case , when S is a positive semidefinite m a trix its corresponding measure is studied in 
Uhlid (Il994l), iDiaz-Garcia and Gutierre3 (Il997t) . iDiaz-Garcia and Gonzale~ Fan'asl (|2005al) 



and IDiaz-Garcia and Gonzalez- Farias! (|2005bf ) under different coordinate systems. 

If A G D^ n then (dA) (the Legesgue measure in ID^J denotes the exterior product of the 
(3m functionally independent variables 



If Hi G V£ _ then 



n 



(d\) = A A dA 



i=l fe=l 



n m 



(HJdHi) = /\ f\ h*dhi. 

where H = (Hi|H 2 ) = (hi, . . . , h m |h m+ i, . . . , h„) S il^(m). It can be proved that this 
differential form does not depend on the choice of the matrix H2 and that it is invariant 
under the transformations 



Hi — > QHP 1; Q G it^(n) and P G it^(m). 



(1) 



When m = 1; Vf „ defines the unit sphere in J™. This is, of course, an (n ~ 1)(3- dimensional 
surface in J™- When m = n and denoting Hi by H, (H*dH) is termed the Haar measure 
on iP(m) and defines an invariant differential form of a unique measure ^ on iP(m) given 

by 



i/(£Dt) 



(H*dH). 



It is unique in the sense that any other invariant measure on il^(m) is a finite multiple of v 
and invariant because is invariant under left and right translations, that is 

i/(QTl) = v{WQ) = v{Wl), VQeil^(m). 

The surface area or volume of the Stiefel manifold Vf n n is 



Vol(Vf n ,J = / (HldHx) = 



Ti[n(3/2] ' 



(2) 



■5 



where r^Ja] denotes the multivariate Gamma function for the space S^(to), and is defined 

by 



r /3 



etr{-A}|A| Q 



-(m-l)/9/2- 



\dA) 



m 

m(ro-l)/3/4 TT pj 



where etr( 



Gross and Richards! (|l987h . 



expftr(-)) 
ardsl ( 

1p 



|o-(i-l)/3/2], 
denotes the determinant and Re(a) > (m — l)/3/2, see 



Dimitriul (J2002, p. 35) first assigned the responsibility of the calculating some Jacobianos 



to researchers of random matrix theory, although in the end this fact is not especially 
important, as the Jacobian of the spectr al facto r isatio n (and of SVD among others) for the 
symmetric matrix was first computed by IJames (Il954h and not by Wigner (1958). 

There exist an extensive bibliography on the computation of Jacobians of matrix trans- 
formations and decompositions in both random matrix theory and statistics. We now sum- 
marised diverse Jacobians in terms of the parameter /3, some based on the work of Dimitriu 
Dimitriu (2002), w hile ot h er results are prop o sed a s extensions of r eal or complex c a ses, see 
Deemer and Olkinl (Il95ll).l01kinl (Il953l). Rovl (fl957h . I James] (fl95l). I Anderson] (ll958l).|Jame£ 



is 



1964) iKhatril (Il965l). iMuirheadl (I1982I) . iMethal (|l99ll ). lOlkinl (|2002f ). iRatnaraiah et al 



2005a) and R atnaraiah et all (|2005bf ). An excellent reference for real and complex cases 
Methail (jl997l ). which includes many of the Jac obian that have b een published. Some 



Jacobians in the quaternionic case are obtained in iLi and Xue (|2009h . We also included a 
parameter count (or number of functionally independent variables, #fiv), that is, if A is 
factorised as A = BC, t h en th e parameter count is written as #fiv in A = [#fiv in B] + 
[#fiv in C], see|Dinfiriiu| (|2002j). 



Lemma 2.1. 

and C G £!L , 



Let X and Y e £.^ n , and let Y = AXB + C, where A 6 C^ n , B G C^ m 
are constant matrices. Then 



Lemma 2.2. Let X and Y 6 

constant matrices. Then 



(dY) = |A*A| /3m / 2 |B*B| /3 "/ 2 (dX). 

<Pf„, and let Y = AXA* + C, where A and C G C^ m 



(3) 



(dY) = |A*A| /3 (" l - 1 ^ 2+1 (dX) 



Lemma 2.3. Let J G X^(m) and write J 
G G %y(m), with gn = 1 for all i = 1, . . . ,n 
• parameter count: (3m(m + l)/2 = [(3m 

rn 

(dj) = n 



(4) 

BG where B = diag(6i, . . . , b m ) G £>^, and 
Then 

\- [[3m(m — l)/2] and 



(dB)(dG), 



(5) 



i=l 



whe 



(dG) 



Hj ■ 



A A M 

i<j k—1 

Lemma 2.4 (LU decomposition, Doolittle's version). Let X G Cf^ „ 

. .to and T G X^-(to). Then 
l)/2] + [f3m(m + l)/2] and 



where A G % L (m), with 8u 
• parameter count: ftm 2 



--l,i = 1, 
= [f3m(m 

m 

(dx)=ni 



|/3(m-i) 



(dT)(dA), 



id write X = AT, 



(G) 



G 



where(dA) = AT> 3 ALidS^. 

Lemma 2.5 (LU decomposition, Crout's version). Let X G CJ^ nm , such that X = AT, 

where A £ T^(m) and T G T^-(rn) si/c/i that va — 1 /or all i = 1, . . . , m. Then 
• parameter count: (3m 2 = \(3m(m + l)/2] + [/3m(m — l)/2] and 

m 

m-ni^i^"" ™^)- (7) 



^ere(dT)=A^ALi^- 

As a consequence of Lemmas 12.41 and 12.31 there follows 

Lemma 2.6 ( LDM decomposition). Let X G £? n rn , such that X = Vt/IIH. where \& G 
T^(m) iot£/i -0m = 1 /or a/i i = 1, . . . , m, II = diag(7Ti, . . . , 7r m ) G and H G T^-(m) ™£/i 
— 1 / or all i = 1, . . . ,m. Then 

• parameter count: j3m 2 — [/3m,(m — l)/2] + [[3m] + [(3m,(m — l)/2] and 

m 

(dx) = n \* ll \ 2 ^ m - l \d*)(<m)(d'z). (8) 

Lemma 2.7 (Qi? decomposition). Lei X G i/iere i/ien exisis an Hi G V/^ n and a 



T G T^(m) reaZ ij, > ; i = 1, 2, . . . , q such that X = ELT. Then 

• parameter count: (3mn = [(3mn — (3m{m — l)/2 — m] + [f3m{m — l)/2 + m] and 

m 

(dX) = Y[tl (n - t+1) - 1 (U* 1 dU 1 )(dT). (9) 

i=l 

Now, from Lemmas 12.71 and 12.41 it follows that 

Lemma 2.8 (Modified QR decomposition (QDR)). Let X G C^ n n , then there exist Hi G 
V^„, a diagonal matrix N = diag(ni, . . . , n m ) G J)^, with n\ > ■■■ > n m > and 
f2 G 2^(to) with uja = 1, i = 1, 2, . . . , m suc/i i/iai X = HiNfJ. TTien 

• parameter count: (3mn = [(3mn — (3m{m — l)/2 — to] + [m] + [[3m(m — l)/2] and 

m 

(dX) = 2- m ]Jnf ( " +m " 2!+1) " 1 (H*dH 1 )(dN)(dO). (fO) 

Lemma 2.9 (Polar decomposition). Lei X G n and write X = PiR, with Pi G n , 
and R G Then 

• parameter count: (3mn = [(3mn — (3m{m — l)/2 — to] + [f3m(m — l)/2 + m] and 

m 

(dX) = |R|/»(»-m+l)-l JJ( di + d J ) /5 (dR)(PtdPi), (f 1) 

where R = QDQ* is i/ie spectral decomposition o/R, Q G il^(m) and D = diag(di , . . . , d m ) G 
with di > • • • > d,„ > 0. 

Lemma 2.10 (Singular value decomposition, SVD). Let X G C@ n n , such that X = ViDW* 
with Vi£V^„, We ttP(m) and D = diag(di, • • • , d m ) G D^, di > • • • > d m > 0. Then 
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• parameter count: (3mn — [(3mn — f3m(m— 1)/2 — m — (J3— l)m\ + [m] + [[3m(m— l)/2+m] 
and 

m m 

(dX) = 2- m TT T Y[dP {n - ,n+1) - 1 Y[(d* - d 2 j f(dB)(V* 1 dV 1 )(W*dW), (12) 

i— 1 i<j 

where 



T = < 



o, 


/? = 


—m, 


/? = 


-2m, 


/? = 


—4m, 


/? = 



Lemma 2.11 (Cholesky's decomposition). Let S £ ^P^> and write S = T*T ; where T £ 



lfj(m) with tu > 0, i = 1, 2, . . . , m. Then 



• parameter count: (3m[m — l)/2 + m = [3m{m — l)/2 + m and 

m 

(dS) = 2 m Y[4 {m - t)+1 (dT). (13) 

i=l 

From Lemmas 12.1 II and 12. 41 it follows that 

Lemma 2.12 ( L*DL decomposition). Let S £ and write 5 = Cl*Ofl, where 17 £ 
T^(m) with uja = 1, i — 1,2, ... , m and a diagonal matrix O = diag(oi, . . . , o m ) £ 33^' 
ww'f/i oi > ■ • • > o m > 0. Then 

• parameter count: (3m(m — l)/2 + m = [/3m(m — l)/2] + [to] and 

m 

(ds)=n«f (m ^w(rfo). (u) 

Lemma 2.13 (Hermitian positive definite square root). Let S and R £ suc/i i/iai 
S = R 2 . T/ien 

• parameter count: (3m(m — l)/2 + to — [3m(m — l)/2 + m and 

<? 

(dS) = 2 m |R| J](d t + d J ) /3 (dR), (15) 

where R = QDQ* is i/ie spectral decomposition o/R, Q £ il^(m) and D = diag(di, . . . , d m ) £ 
35*, wii/i d\> ■■■ > d m > 0. 

Lemma 2.14 ( Spectral decomposition). Let S £ T/ien </ie spectral decomposition 

can be written as S = WAW*, where W £ it^(m) and A = diag(Ai, . . . , A m ) £ 2)^, wit/i 
Ai > • • • > A m > 0. Then 

• parameter count: (3m[m — l)/2 + m = [/3to(to + l)/2 — m — (/3 — 1)to] + [to] and 

m 

(dS) = 2-" 1 tt t YliX, - A 3 ) /3 (dA)(W*dW), (16) 

i<j 

where t is given in Lemma \2.1(h 

Finally, by combining Lemmas HHUl O and [Ml with Lemmas I2TT1 I2TT51 [2TTT1 I2TT21 

respectively, the following result is obtained. 
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Lemma 2.15. Let X e C^ n , and write X = ViDW* (SVD), X = ViR (Polar decom- 
position), X = ViT (QR decomposition) or X = ViNfi (modified QR decomposition) and 
let S = X*X e q}^. Tften 

(dX) = 2- m |S| /3 ("- m+1 )/ 2 - 1 (dS)(V*dV 1 ). (17) 

Note that the Lebesgue measure (dS) in jl7]) , can be factorised in terms of Cholesky, L*DL, 
Hermitian positive definite square root or spectral decompositions, thus obtaining this way 
alternative explicit expressions of(dS), give n by H3\), (14 \), (15\) and if 16\) , respectively. For 
the corresponding coordinates and bases, see \Eatori kOOl , Chapt er ol. 



3 Elliptical ensemble 

In this section we employ the nomenclature recently proposed in random matrix theory for 
some of the most commonly-studied random matrices. These terms are Hermite, Laguerre, 
Jac obi and Fourier instead of G a ussian, Wishart, MANOVA ( or beta type I) and c i rcular , 



SCO 



Dimitriu and Edelmanl (|2002h . lEdelman and Raol (|2005h . lEdelman and Suttonl (|2008h 



These alternative names are proposed in view of the fact that in random matrix theory we 
are interested in matrices with joint eigenvalue density proportional to JXiLi w (^i)\A{^)\^ 
where |A(A)| — Yii<j(^i ~ ^j) 1S the absolute value of the Vandermonde determinant and 
w(X) is a weight function of a system of orthogonal polynomials. In particular, when the 
joint eigenvalue density is calculated under Gaussian, Wishart and MANOVA random ma- 
trices, these weight functions correspond to the Hermite, Laguerre and Jacobi orthogonal 
polynomials, respectively. 

Under generalised multivariate analysis, f our classes of ma t rix va riate elliptical distri- 
butions have been defined and studied, see iFang and Zhanel (|l990h . A random matrix 
X G C) n n is said to have a matrix variate left-elliptical distribution, the largest of the four 
classes of matrix variate elliptical distributions, if its density function is given by 

C(m,n) , fo(S" 1/2 (X ~ ^) T © _1 (X - M )£- 1/2 



| S |n/2| |m/2' 



where h is a real function, c(m, n) denotes the normalisation constant which might be a 
function of another parameter implicit in the function ft, £ € £ ^ and fi e $l nXTn . 
This fact is denoted as 

X - ££<S mx „(/u, S, 0, ft,). 

On e might wish to consider situ ations where X has a density function of the following form, 
see 



Fang and Zhanel (jl99dh and lFang and Lil (|l99&h . 



c(TO ' n) ^(s^cx-Mfe^cx-i*)). (18) 



| S |n/2|0|m/2 

This condition is equivalent to considering the function ft as a symmetric function, i.e. 
g\g(AB) — g(BA) for any symmetric matrices A and B. This condition is equivalent to 
that in which ft(A) depends on A only through its eigenvalues, in which case the function 
ft(A) can be expressed as ft(A(A)), where A(A) = diag(Ai, . . . , A m ) and Ai, . . . , A m are the 
eigenvalues of A. Two subclasses of matrix variate elliptical distributions are of particu- 
lar interest: matrix variate vector- spherical and spherical elliptical distributions. For these 
distributions, A(A) = tr(A) and A(A) represent any function of eigenvalues of A, and are 
denoted as X ~ £VS mxn (fi,Y^,®, h) and X ~ £SS mxn (fi, S, 0, ft), respectively. Note 
that matrix variate vector-spherical elliptical distributions are a subclass of matrix variate 
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spherical elliptical distributions. Many well-known distributions are examples of these sub- 
classes; one such is the matrix variate Hermite distribution. Other variants include matrix 
variate vector-spherical elliptical distributions, e.g. Pearson type II, P e arson type VII, Kotz 
type , Bessel and Logistic, among many others, see iGupta and Vargal (|l993l ). Further cases 
are those of matrix variate spherical elliptical dis tributions, e.g. Pea rson Type II, Pearson 
type VII and Kotz type, among many others, see iFang and Li dl999h . 

When in matrix variate vector-spherical and spherical elliptical distributions it is as- 
sumed that S = I m , = I n and fi = nxm we obtain the matrix variate vector-spherical 
(denoted as V5 mxn (0, I m , I„, h) = VS mxn (h)) and spherical distributions (denoted as 
(0,I m ,I„, h) = SS mxn (h)), which are known to be invariant distributions under or- 
thogonal transformations, because X and QXP have the same distribution when P 6 il 1 (m) 
and Q € it^n). In particular, note that if ir is a random variable with vector-spherical or 
spherical distribution it is denoted as x ~ VS(0, 1, h) = VS(h) or x ~ SS(0, 1, h) = SS(h), 
respectively. 

We now define spherical and generalised Laguerre ensembles based on vector-spherical 
and spherical distributions. As can be seen, these ensembles contain as particular cases the 
classical Hermite, Laguerre, Jacobi and Fourier ensembles, as well as many o thers that have 
been studied in the literature of random matrix theory, see iForrester (200J). 



3.1 Vector-spherical and spherical random matrices 

V (n,m, h) and SS@ (n, to, h) are nxm matrices of non-correlated and identically dis- 
tributed (n-c.i.d.) with entries V<S^(0, 1, h) and 55^(0, 1, h), respectively. 



1. If A G Ct 



is an n x m 



vector-spherical random matrix VS^(n,m,h) then its joint 



element density is given by 



c /3 (m,n)/i(tr(A*A)) 



(19) 



where 



2. If A e£^ 



T[/3mn/2 
c 1 (m,n) — ■= 777 




(3mn—l 



h(u 2 )dv 



is an nxm 



spherical random matrix SS@(n, m, h) then its joint element 



density is given by 
where dP(m 7 n) is given by 
rgJ/3n/2]rgJ/3m/2 

_m(n+l)/3/2+T 



d l3 (m,n)h(X(A*A)) 




■J |L|/9(«-^-+i)/2-i|A(L)|^/ i (L)(dL) 



(20) 



(21) 



h>—>l m >0 



with L = diag(?i, . . . , l m ) and U is the ith eigenvalu e of (A* A) and r is giv en in 
Lemm a 12.101 for the re al case see IFang et al. (|l990l) . IFang and Zhanl (l990) and 



Gupta and Vargal (|l993l) 



Note that A and QAP have the same distribution when P £ it^(m) and Q G it^(n). 
That is, the distributions of VS^(n, m, h) and SS@(n, m, h) are invariant under orthogonal 
((3 = 1), unitary ((3 = 2), symplectic (/? = 4) and exceptional type (f3 — 8) transformations. 

Remark 3.1. Note that the unique case in which the elements ay of the random matrix 
A are i.i.d. is when V5^(0, l,h) = SS p (0, 1, h) = AfP(0, 1), this is when q„- has a standard 
normal (Gaussian or of He r mite) distribution, see lFang and Zhan g (1990, Theorem 2.7.1, p. 
72) and lGupta and Vargal (|l99a Theorems 6.2.1 and 6.2.2 p.193). 
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Some particular vector-spherical random matrices and their corresponding density func- 
tions are summarised in the following result, where u = tr A* A). 

Corollary 3.1. Let A 6 C^ n n with density function U9)) then 

1. (random matrix of Hermite) If h(u) = exp(—(3u/2), it is said that A is a random 
matrix of Hermite with density function 

J etr{-/3A*A/2}. 



(2Trf3- 1 ) l3mn/2 

2. (T type I random matrix) If h(u) = (1 + u/v)-^ mn+ ^ 2 \ it is said that A is aT 
type I random matrix with density function 

where v > 0. This distribution i s also term ed the Pearson type VII matrix v ariate 
distribution, see Guvta and Farad tl99& ) and \ Diaz- Garcia and Gutierrez 1 200& ). 

3. (Gegenbauer type I random matrix) If h(u) = (1 — u)^ l3q , it is said that A is a 
Gegenbauer type I random matrix with density function 

r^mn/2 + /? g + l] . A * A1 - gq 

T T 0m n /2 TP[(3q+l] ^ HA A; , 

where q > —1 and tr A* A < 1. This distribution is known in stati stical bibliography 
as the i nverte d T or Pear son type II matrix variate d i stribu tion, see Guvta and Varaci 
' 199j ), \Presl Xl98& ) and lDiaz- Garcia and Gutierrek \200A ). 



For the spherical random matrices we have. 
Corollary 3.2. Let A 6 Cf^ n n with density function i!9\) then 

1. (T type II random matrix) If h(\(A* A)) = |I m + A*A\-> 3( - n+ ^/ 2 , it is said that 
A is aT type II random matrix with density function 

r^[/3(n + v)/2] + A * A |-0(n+!/)/2 



(i)^"/2 r ^[^/2] 



where v > m. This distri bution i s also t ermed the T o r Pearson type VII ma- 
tricvariate distribution, se e Dickeil 1 196$ ), \Prest (198& ), \Fana and L\ (199& ) and 



\Piaz- Garcia and Gutierrei 1 200, 
2. (Gegenbauer type II random matrix) If 

h(X(A*A)) = \I- A*Af^- m+1 ^ 2 -\ 
it is said that A is a Gegenbauer type II random matrix with density function 
r&[/3(n + i/)/2] 



7r/W 2 rf„[/^/2] 



|j _ A * A |0(«/-m+l)/2-l ) 



where v > (m — l)/2 and A* A 6 This distribution is known in stat istical 

bibliog r aphy as the inverted T or Pearson type II matricvari a te dis tribution, see Press] 
■ l9S.l,.\Fana and Li \l999\ ) and \Diaz- Garcia and Gutierrek \200& ). 
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3.2 Construction of the vector-spherical and spherical random ma- 
trix ensembles 

The vector-spherical and generalised vector-spherical-Laguerre ensembles are constructed 
from VS 13 ^, m, h) as follows. 

Vector-spherical orthogonal ensemble (VSOE): symmetric mxm matrix obtained 
as (A + A T )/2, where A is FS' 1 (m, to, h). The diagonal entries are n-c.i.d. with distribu- 
tion V«S (0, l,h), and the off-diagonal entries are n-c.i.d. (subject to the symmetry) with 
distribution VS^O, 1/2, h). 

Vector-spherical unitary ensemble (VSUE): Hermitian mxm matrix obtained as 
(A + A*)/2, where A is VS 2 (m, to, h) and * denotes the Hermitian transpose of a complex 
matrix. The diagonal entries are n-c.i.d. with distribution V«S (0, l,h), and the off-diagonal 
entries are n-c.i.d. (subject to the symmetry) with distribution VS 2 (0, 1/2, h). 

Vector-spherical symplectic ensemble (VSSE): quaternionic Hermitian (self-dual) 
mxm matrix obtained as (A + A*)/2, where A is VS 4 (m, m, h) and * denotes the quater- 
nionic Hermitian (or dual) transpose of a quaternion matrix. The diagonal entries are 
n-c.i.d. with distribution V«S (0, 1, h), and the off-diagonal entries are n-c.i.d. (subject to 
the symmetry) with distribution V<S 4 (0, 1/2, h). 

Vector-spherical exceptional type ensemble (VSETE): octonionic Hermitian mx 
to matrix obtained as (A + A*)/2, where A is VS (m, m, h) and * denotes the octonionic 
Hermitian transpose of a octonionic matrix. The diagonal entries are n-c.i.d. with distri- 
bution VS (0, l,h), and the off-diagonal entries are n-c.i.d. (subject to the symmetry) with 
distribution VS 8 (0,l/2,/i). 

Analogously, by replacing VS and LS by SS and SS, respectively, in the four previous 
definitions we obtain the Spherical orthogonal ensemble (SSOE), Spherical unitary 
ensemble (SSUE), Spherical symplectic ensemble (SSSE), and Spherical excep- 
tional type ensemble (SSETE). 

Similarly, generalised vector-spherical-Laguerre and spherical-Laguerre ensembles can be 
defined as follows. 

Generalised vector-spherical-Laguerre (V S L 13 (n, m, h), n > m): Symmetric/ Her- 
mitian/ quaternionic Hermitian/ octonionic Hermitian mxm matrix which can be obtained 
as A* A, where A is VS^(n, to, h). 

Again, by replacing VS@ by SS 13 in the previous definition we obtain the Generalised 
spherical-Laguerre (SL /3 (n,m,h), n > m). 

As shown later, the Jacobi ensemble is a particular case of VSL^{n, to, h) or SL^(n, to, h) 
and the Fourier ensemble is a particular case of SS^{m, to, h). 

3.3 Computing the joint element densities 

Let A be an to x to matrix from the vector-spherical ensemble, that is 




(22) 



Then, it is straightforward to see that its joint element density is 



2 ™(™- 1 W 4 c /3(TO,TO)/i(trA 2 ) , 



(23) 



where 



cf(m, to) 



r[(TO + TO(TO- l)/3/2)/2] 
2 7r (m+m(m-l) / 3/2)/2 




',11 



1 
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Remark 3.2. In the general case, from JISJ) we have |£|£V 2 |0|' 3m / 2 = |0 ® S^/ 2 . Under 
(l22|) , ® X is a diagonal matrix with m-times ones and m(m — l)/2-times 1/2, considering 
that A is symmetric, from which |0 <X> E^/ 2 = (i/2)" i ( m - 1 )' 3 / 4 . 



Similarly, let A be an m x m matrix from the spherical ensemble, that is 

SSfi(0,l,h), i = j, 
SS (0,1/2, h), i>j. 



(24) 



And then it follows that its joint element density is 

y»(«- 1 WV(m,m)/i(A(A 2 )), (25) 

where 

d "M = 2m jS/Z 2 ^r { J - J |A(L)|^(L 2 )W| , (26) 

Kh>->l m >0 ) 
with L = diag(^i, . . . , l m ), where h is the ith eigenvalue of A and r is as given in Lemma 

EH 

Consider the generalised vector-spherical-Laguerre ensemble VL(n,m,h) = S = A* A, 
where A = VS(n,m,h). Its joint element den sity can be computed in diverse ways. For 
example, following the approach of iHera (jl955l ), let A = ViR be the polar decomposition 
of A, then S = A*A = R 2 , from Lemma [Tfl and (fl"9)) . the joint density of S and Vi is 

c /3 (m,n)/ l (tr(S))2- m |S|' 3( "- m+1 )/ 2 - 1 ( ( iS)(V^V 1 ). 

On integrating with respect to Vi using @ the marginal density of S is found to be 



^V^ |s p M+ l)/2-l 



h(tv(S)). (27) 



Note that the same result can be obtained from Lemma [T"7l and taking into account the SV, 
QR or MQR decomposition instated of the polar decomposition in the previous procedure. 

As in the generalised vector-spherical-Laguerre ensemble case, it is obtained that the 
joint element density of the generalised spherical-Laguerre ensemble is 

dP(m, n)^ mn ' 2 | S |£(n-m+l)/2-l fa(A(s) y (2g) 
V m [pn/2] 



As was done for Corollarv l3.1i some particular vector-spherical random matrix ensembles 
and their corresponding density functions are summarised in the following result. 

Corollary 3.3. Let A S with density function \23\) then 

1. ( Classical Hermite ensemble) Ifh(u) — exp(— /3u/2), it is said that A is a Hermite 
ensemble with density function 

1 



2>n/2 (• 7r ^-l-) m / 2 + m ( m - 1 )/ 3 / 4 



etr{-/3A 2 /2}. 



2. (T type I ensemble) If h(u) = (1 + u / v )-m/2-p(m(m-\)/2+ v )/2 ^ it is said that A is 
a T type I ensemble with density function 

2m ( m -l)0/4 rm /20 [(3{m{m _ 1)/2 + v)/2] /2 ^ (m[m _ 1)/2+v)/2 
n m/2+m(m-l)f3/4Y0[^/2] { ' 

where v > 0. This ensemble could be termed the Pearson type VII matrix variate 
ensemble. 
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3. (Gegenbauer type I ensemble) If h(u) = (1 — u) ^ q , it is said that A is a Gegen- 
bauer type I ensemble with density function 

2 m(m-i)/3/4 r /3[ m / 2 + m(m - l)/3/4 + /3g + 1] M t _ 2 ^ 



7r m/2+m(m-l)/3/4p/3[ / g 9 + ±] 



-(1 - tr A 2 )" 



where q > — 1 and tr A 2 < 1. In i/iis case, £oo, i/ie ensemble might be termed the 
inverted T or Pearson type II matrix variate ensemble. 

For the spherical random matrix ensembles we have the following. 

Corollary 3.4. Let A 6 mf/i density function i25\) then 

1. (T type II ensemble) If h(X(A 2 )) = \I m + A 2 \-^ n+u ^ 2 , it is said that A is a T 
type II ensemble with density function 

2 m(m- lWV( mjm )| Im+A 2|-/J(n+^/2 j 

where v > m. This ensemble could possibly be termed the Pearson type VII matric- 
variate ensemble. 

2. (Gegenbauer type II ensemble) If h(\(A 2 )) = \I m - AY (l/-m+1)/2-1 , it is said 
that A is a Gegenbauer type II ensemble with density function 

2 ™(™-iWV(m,m)|I m - AY (, '- ro+1)/2 - 1 , 

where v > (m — l)/3 and A* A 6 ^3^. TTiis ensemble could be termed an inverted T or 
Pearson type II matricvariate ensemble. 

Analogously, some particular cases for vector-spherical and spherical-Laguerre random 
matrix ensembles are obtained in the following two results. 

Corollary 3.5. Let S £ ^ n with density function \27^ then 



1. (Classical Laguerre ensemble) If h(u) = cxp(— /3u/2), it is said that S is a La- 
guerre ensemble with density function 

1 , S| /3(n - m+1)/2 - 1 etr{-/3S/2} 



(2/3- 1 f mn/2 r^mn/2] 
where n > (m — l)/3. 

2. (T-Laguerre type I ensemble) If h(u) = (1 + it/zy) - ^"^)/ 2 , it is said that S is a 
T -Laguerre type I ensemble with density function 

Y f3 [f3{mn + v)/2] gi^n-m+n^-i^ + tr g\-/j( n +i/)/2 
TP[/3u/2]Ti[l3n/2Y 

where v > and n > (m — 

5. (Gegenbauer- Laguerre type I ensemble) If h(u) = (1 —u)^ l3q , it is said that S is 
a Gegenbauer- Laguerre type I ensemble with density function 

TP[f3mn/2 + gg + l] | S |/3(n-m+i)/2-i (1 trg) g q 
r/3[/? g + l]lt[/?n/2] ' 
where q > — 1, n > (m — l)/?/2 and tr S < 1. 
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For the spherical-Laguerre random matrix ensembles we have the following. 

Corollary 3.6. Let S G with density function i28\) then 

1. (T-Laguerre type II ensemble) If h(X(S)) = |I m + S\-^ n+u ^/ 2 , it is said that S is 
a T-Laguerre type II ensemble with density function 

1 |g|/3(n-m+l)/2-l|j , g|-^(n+f)/2 



B&\Pm/2,pn/2] 

where v > (to — 1)0 and n > (m — l)/3. XTiis dis t ributio n is also known as the 
Studentised Wishart distribution, see lOlkin and Rubin\ tl96A) . 

2. (Gegenbauer-Laguerre type II ensemble) If h(X(S)) = |I m - } a 

is said that S is a Gegenbauer-Laguerre type II ensemble with density function 

1 'g|/3(n-m+l)/2-l|j m _ g |/3(i/--m+l)/2-l 



where v > (m — and n > (to — l)/3. 
Finally, note that B^[a, b], defined as 



is tfte multi variate beta function, where Re(a) > (m — l)/3/2 and Rc(6) > (m — l)/3/2, see 

S3 



Remark 3.3. 1. Note that, f rom flj, Hi G V^„ is a spherical random matrix. More- 
over, (see Fang and Zhand ( 1990l Lemma 3.1.3(iii), p. 94)), the differential form of its 
density function is 

Vol (yi,n) 2™7T^"/ 2 

In the context of random matrix theory, Hi is a Fourier random matrix. Now, let 
itf (m) be the group of orthogonal Hermitian matrices H. Then H is a Fourier ensemble 
with a differential form of its density function 

1 (H.H) = ^»&(H^). (29) 



Vol (m)) 2™ 7 rW + r (rW 2 + 1]) 

The value of Vol (to)^ is found in next subsection. 

2. Some g(A) functions, where A is a spherical random matrix, have invariant distribu- 
tions under the corresponding class of spherical distribution, under certain conditions; 
in other words, Y = g(A) has the same distribution for each particular spherical dis- 
trib ution. An analog o us sit uation is true for the class of vec tor-spherical distribution 



see 



Fang and Z hang I (|1990L Section 5.1, pp. 154-156) and iGupta and Vargal (|1993l . 



Section 5.1, pp. 182-189). In particular note the following: let A G „, be an n x to 
spherical random matrix SS^(n,m, h), such that A = (Af|A^) T with A, ; G 
Ui > m i = 1, 2 and n = rii + n 2 . Then 
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(a) the T type II random matrix defined as 



Ai(A£A 2 ) 



-1/2 



(b) the Gegenbauer type II random matrix defined as 

R = A 1 (A 1 A 1 +A;A 2 )- 1 / 2 , 

(c) the classical Jacobi random matrix defined as 

B = (A*Ai + A^A 2 )- 1 / 2 (A 1 A 1 )(AtA 1 + A^A,)" 1 / 2 , 

(d) and the F random matrix defined as 

F = (A*A 2 )- 1 /2(AtA 1 )(A*A 2 )- 1 /2, 



are invariant under the class of spherical distributions, see IFang and Zhang (ll990L 
Section 3.5, pp. 110-116). Analogous results are concluded for the class of vector- 
spherical distributions, see lGupta and Vargal (|l993l Theorem 5.3.1, p. 182). Moreover, 
the classical Jacobi and F random matrices have the same j oint element distr i bution 
if A is a vector-spherical or spherical random matrix , see IFang and Zhang (|1990L 
Theorems 3.5.1 and 3.5.5) and Gupta and Varga ( 19931 Theorem 5.3.1). 



3. Observe that the Gegenbauer-Laguerre type II ensemble is indeed the classical Jacobi 
random matrix ensemble. Also note that, if F is an F random matrix, then B = 
I m — (I m + F) _1 , and if B is a classical Jacobi random matrix ensemble then F = 
(I m — B) _1 — I m . Therefore F can be termed a modified Jacobi random matrix 
ensemble. Similarly, the T-Laguerre type II random matrix ensemble is indeed the 
modified Jacobi random matrix. 

4. There remain two final remarks: 



(a) As lEdelman and Suttonl (|2008f) observed, the classical Jacobi random matrix can 
be obtained from the Fourier random matrix. Moreover, this procedure can be 
applied to any function g(A) 7 invariant under the class of spherical distri butions. 
Form the statistician's point of view, this approach has been used by IKhatri 
(|l970h and ICaded (|l996 h . in the real case. The latter studied the T, inverted T, 
F and beta distributions. 

(b) Finally, note that the matrix factorisation associated with classical and modi- 
fied Jacobi ensembles could be the singular value dec ompositions, but applied t o 
Gegenbauer and T random matrices, respectively, see lEdelman and Rao (|2005h . 



3.4 Joint eigenvalue densities 

Theorem 3.1. Let A be an m x m matrix from the vector-spherical ensemble. Then its 
joint eigenvalue density function is 



rP(m T1 N lO"i(m-l)/3/2_/3m 2 /2+r m ( m \ 



r&[/3m/2] 



(30) 



*<3 



where A = diag(Ai, . . . , A m ). 



16 



Proof. From (|23[) and Lcmma r2.141 the joint density of W 6 UP{m) and A = diag(Ai, . . . , A m ) 
with A = WAW* is 

m 

2 ™('"- 1 )^/V(m,m)/i(trA 2 )2- Tn 7r T JJ(A i - Xjf(dA)(W*dW). 

i<j 

The marginal density desired is found by integrating over W S ll' 3 (to) using @. □ 
Similarly, 

Theorem 3.2. Let A be an m x m matrix from the spherical ensemble. Then its joint 
eigenvalue density function is 

/ l (diag(A?,... I A^ l )) (31) 

where A = diag(Ai, . . . , A m ). 

Proof. The proof is analogous to that given for ([3H)1 ; note however that /i(A(A 2 )) = 
MA((WAW«) 2 )) =/i(A(A 2 )) =/ l (diag(Af,... ) A^)). □ 

Theorem 3.3. Lei S be an m x m matrix from the vector-spherical-Laguerre ensemble. 
Then its joint eigenvalue density function is 

rn 



rP(m «W 3m (™+ m )/ 2 + r m m 
?m - iiW TT^(n-m+l)/2-lTT( A . _ Aj .)/3 

Mf>], (32) 



r^[/3n/2]r^[^m/2] fj 



vi=l 



where A = diag(Ai, . . . , A m ). 

Lroo/. Let S = WAW*, from ([27]) and Lemma l2~T4l then the joint density of W e it (to) 
and A = diag(Ai, . . . , A m ) is 

° [m ^ n)7T |S|' 3 ("- m+1 )/ 2 - 1 /i(tr(S))2- m ^ JJ(Ai - A J ) /3 (rfA)(W*dW). 

The marginal density desired is found by integrating over W S Ur(m) using @. □ 
Similarly, 

Theorem 3.4. Let S be an m x m matrix from the spherical- Laguerre ensemble. Then its 
joint eigenvalue density function is 

= c (m,n)7r -pr ^(„_ m+ i)/2-i TT (Ai _ x ) 

Ti[(3n/2}ri[f3m/2] l\ j£ 

x /i(diag(Ai,...,A ro )), (33) 

where A = diag(Ai, . . . , A m ). 
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Given the differential form of the joint element density function of the Fourier ensemble, 
the corresponding joint eigenvalue density function cann ot be obtaine d by applying the 
Theorem 13.21 This densi ty function has been obtained by iMetha ( 1991 . Lemma 10.4.4, p. 
198) and lForresteil (|2009l Proposition 2.3, p. 61). Based on the wedge product, this is now 

obtained and Vol ^]3g(m)J is calculated indirectly, see (|29p. 

Let H e ?pf(m), then there exist U G ^(m) such that H = UEU*, wher e E = 
diagfe xpfi fli ), . . . exp ( ifl m )) , exp(i$i) are m complex numbers on the unit circle, see IMetha 
(1991) and iForresterl ( 2009f ). Then the joint density function of 9\, . . . , 6 m is 

c^(m) JJletpW - exp(i0j)f , #j e (-tt.tt). 

l<3 



From ([29 



dF H (H) 



(HdH) 



Vol (?pf(m)) 

Let H = UEU*, then by IForresterl (|2009l . Exercise 2 (iii), p. 63), recalling that H = H* 
(HdH) = J]|exp(^)-exp(%)l /3 ( A ') ( U * dU ) 

i<7 \i=l / 

Note that it is necessary to divide the measure by 2" m 7r T to normalise the arbitrary phases 
of the m elements in the first row of U, where r is given in Lemma 12.101 Then the joint 
density of diag(#i, . . . 6 m ) and U is 



Vol 



(?pf(m) 



U I exp(i0,) - exp(^)|^ (f\6i) (U*dU). 



i<j \z=i 
By integrating with respect to U we obtain the marginal of B\ , 
2-' m n T Vol (^(m)) 



Vol 



(?pg(m)) 



niexp(^)-exp(^)|^/\i 



(34) 



z=i 



From the Morris integral Forresterl ( 20091 Equation (3.4), p. 122) 

/7T /*7T m 

■■■/ IJ|exp(t0 J )-exp(t0 3 -)| /, /\0, 

(27r)- m r[/3m/2 + 1] 



/=i 



then 

However, from (|34[) we have 
c"(m) = 

Therefore 



(r[/3/2 + i]) m ' 

r[/3m/2+ 1] 



cP(m) 



(r[/3/2 + 1])- ■ 

r[^m/2 + 1] 2-" l ir T Vol (5p^(m)) 



Vol 



(r[/?/2 + !])"> Vol(^(m)) 
2" l 7r' 3m2 / 2 + T (r[/3/2 + 1])" 



(^f(m) 
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Table 2: Kernels of the joint eigenvalue density function of all the particular ensembles 



Ensemble 

\ t e(a.b) 


Kernel 


Hermite 

Xi E( — oo, oo) 


2 * — ' 

I, i=l ) i<j 

/ m \ -/3[m(m-l)/2+H/2-m/2 m 

V i=l / «j 
/ m \ ^9 m 

m m 


T type I 

A^ E( — 00,00) 

Gegenbauer type I 

Ai6(-l,l) 


T type II 

A^ G(— 00,00) 


11(1 + IJ (A4 _ Ai) /» 

i—1 z<j 
m m 


Gegenbauer type II 

A*6(-l,l) 


JJ (1 _ A 2 ) ^- m+1) /2-l-Q (A ._ A . )/3 

i—1 z<j 


Laguerre 

A;G(0,oo) 

Modified Jacobi type I 

A;e(0,Oo) 

Jacobi type I 

A,e(o,i) 


m ( q m ^ m 
Z — ' 

8=1 k i=l ) i<j 

m / m \ -/3("+f)/2 m 

TT A /5(«-™+i)/2-i [ uVi. TTfAi - A„y 

ll^i I Z_— / I 11^ 4 J ' 

i=l \ i—1 J i<j 

i=l \ i—1 / i<j 
m m 


Modified Jacobi 

A;G(0,oo) 


J] A f ( "- m+1)/2 - 1 (l + A,;)-' 3 ^/ 2 JJ(Ai - Xjf 

i—1 i<j 
m m 


Jacobi 

A;G(0,1) 


J] A f («-™+D/2-l (1 _ Ai) /3(,-m+l)/2-l JJ (A . _ Xj)P 

i—1 i<j 
m 


Fourier 

Aie(-7T,7r) 


JJ|exp(^)-exp(^)| /3 



Table [2] summarises the kernels of the joint eigenvalue density function of all the par- 
ticular ensembles studied in the subsection 13.31 using the names proposed in Remark 13.31 
The corresponding normalisation constants are obtained from their joint element density 
functions and from Theorems 13-11 13. 2i 13.31 and 13.41 

Finally, another interesting property of the vector-spherical random matrix ensemble is 
that its joint normalised eigenvalue density function is invariant under all classes of vector- 
spherical distribution, where the normalised eigenvalue can be defined as <5; = XJr. In 
particular, r can be defined as 

/ m \ 

r =II> 2 ] ■ 



This problem has been st udied in the contex t of sh ape theory and the corresponding density 
function was obtained bv iDfaz-Garcia et~al. ( 20031 ). 
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4 Multivariate statistics 

Many areas of multivariate statistics (and of statistics in general), such as multiple time 
series or e conometr i cs, ca n be enriched with the use of the tools and ideas of random matrix 
theory, see lHannan ( 197dh . However, in this section these tools and ideas are used to develop 



a unified theory of multivariate distributions for normed division algebras. 

Theorem 4.1. Let Z G C^ nn be a matrix variate left elliptical distribution with density 
function 

c /3 (m,n)h(Z*Z). 

Therefore if X = AZB + fi, with A g £J n , Be £^ m and fi G C^ n n , constant matrices 
such that A* A = G qjg and B*B = £ = (S 1 / 2 ) 2 g fAen 



c^(m, n) 



|S|/3n/2|0|/3m/2 

where c^(m,n) is a normalisation constant. 

Proof. The proof is analogous to that given in the real case, but considering the Jacobian 
of the transformation Y = AXB + /j, defined by Lemma \2. II □ 



Corollary 4.1. Under the hypothesis of Theorem \4-l\ 

1. for the matrix variate spherical elliptical distribution, its density function is 

—^f^£L—h (a (5r x (x - i ti)*0- 1 (x - ti))) , 

this fact being denoted as X ~ SS£^ xm (fi, S, ©, h). 

2. for the matrix variate vector-spherical elliptical distribution, its density function is 

—^h^L—h (tr (S- X (X - jl i)*0- 1 (X - fi))) , 

this fact being denoted as X ~ ViS£^ xm (/i., S, 0, h). 

3. and for the matrix variate normal distribution, its density function is (taking h(u) 
exp(-/?u/2)J 



1 

etr 



(2ltl3- X f mnl ' 1 |£|/W2|0|/3m/2 

denoting this fact osX~ A/jf xm (/x, S, ©) 



j-^s-^x-^^-^x-^l 



Many other particular vector-spherical or spherical elliptical distributions can be ob- 
tained by simply specifying the function h(-) in a similar way to Corollaries 13.31 and 13.41 

Elliptical and, in particular, normal symmetric random matrices have received less at- 
tention in multivariate statistics, but analogous results can be obtained in a similar way 
using Lemma [2~2l 



Theorem 4.2. Let us define S = X*©"^ G the 
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1. i/X ~ SS£ „ xm (0, S, 0, ft), the density function ofS is 

dP{m, fl)^™' 2 |g|ff(n-m+l)/2-l fe(A(S -l S))- 

r^[/3n/2]|S|^/2' 

Tftis distribution is known as the spherical- generalised- Wishart distribution and it is 
denoted as S ~ «S£W^(n, S, ft). 

2. i/X ~ ViSf „ xm (0, S, ©, ft), ifte density function of S is 

d0(m,7i)7r /3mn / 2 



r^ l [/3n/2]|S|' 3 «/2 



|S| /3 ("- ro+1 )/ 2 - 1 /i(tr(S- 1 S)). 



Tftis distribution is known as the vector-spherical-generalised-Wishart distribution and 
it is denoted as S - V<S5W^(n, S, ft). 

Proof. Let S = X*© _1 X, and note that 

Y = 0-V2 X ~ 55^^(0, S,I„, ft) 

with (0- 1 / 2 ) 2 = 0. Therefore S = X*0 *X = Y*Y. The desired results are follow from 
and (|2"T|) . respectively. □ 



Corollary 4.2. Under the hypothesis of Theorem \4-S\ it follows that, 
1. the Wishart random matrix has a density function 

1 



|g|/J(n-m+l)/2-l etr {_ / 3 S -l S /2} 



(2/3^ 1 ) /3 " m/2 rf„ [0mn/2] | £ | /W 2 
where n > (m — and iftis is denoted as S ~ W^(n, £).; 
^. £/ie matrix variate beta type I distribution has a density function 

' -|S| /3 ( n_m+1 ^ 2_1 |I m - s| /3 ( iy - m+1 )/ 2 ~ 1 



B^[/3n/2,/3i//2]' 

where v > (m — l)/3 and rt > (m — and this is denoted as S ~ £>Z^(n, v); 
3. and the matrix variate beta type II distribution has a density function 

1 -isi^"-" 1 ^/ 2 - 1 !!™ + s\-^ n+v ^ 2 , 



Bi[pm/2, f3n/2] 

where v > (m — l)/3 and n > (to — and this is denoted as S ~ BT2^ n (n, v). 

Proof. The desired results are obtained from Theorem l4.2l directlv. after defining ft(tr(S _1 S)) 
etr{-/3S- 1 S/2}, ft(A(S" 1 S)) = |I m - s\P(v-m+i)/2-i (s = ^ and fe ( A ( S -i S )) = | Im + 
s \-0(n+v)/2 ( S = reS p ect i V ely. □ 
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Conclusions 



Although most results about Jacobians are known in the context of random matrix theory, 
in general for statisticians they are less familiar, as are the corresponding technical tools 
in the context of normed division algebras. Thus the importance of addressing this area of 
study. 

In the field of random matrix theo ry, various classe s of ensembles are proposed, including 
many of those studied previously, see Forrester ( 20091 . Section 4.1.4, p. 177). What is most 
important is that these new classes of ensembles contain many other ensembles of potential 
interest, which may enable us to study phenomena and experiments under more general 
conditions. 

Analogously to current random matrix theory, see lEdelman and Raol (|2005h . in the 
present study we propose a unified theory of matrix variate distribution for normed division 
algebra, that is, for real, complex, quaternion, and octonion cases. 
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